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2 Abstract 

^ The quantum behavior of superconducting qubits coupled to resonators is very similar to that of 

55 atoms in optical cavities [1, 2], in which the resonant cavity confines photons and promotes strong 

light-matter interactions. The cavity end- mirrors determine the performance of the coupled system, 

s . . . 

I with higher mirror reflectivity yielding better quantum coherence, but higher mirror transparency 

I ... _e. .„.o,, ... . . . . ... 

I mirror transparency, enabling switching between long photon lifetime during quantum interactions 

T— H and large signal strength when performing measurements. Here we demonstrate the superconduct- 

> 

ing analogue, using a quantum system comprising a resonator and a qubit, with variable coupling 
^ to a measurement transmission line. The coupling can be adjusted through zero to a photon emis- 

00 sion rate 1,000 times the intrinsic photon decay rate. We use this system to control photons in 

o 

CN coherent states as well as in non-classical Fock states, and dynamically shape the waveform of 

released photons. This has direct applications to circuit quantum electrodynamics [3], and may 
enable high-fidelity quantum state transfer between distant qubits, for which precisely-controlled 

H 

waveform shaping is a critical and non-trivial requirement [4, 5]. 
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Superconducting resonators play a central role in quantum information technology. Ap- 
plications include the synthetic generation and storage of photon quantum states [6-8], 
quantum memories for quantum computation [9] , and dispersive measurements of supercon- 
ducting qubits [10, 11] as well as defects in diamond [12, 13]. Resonators with low internal 
losses are typically desirable, but the resonator's coupling strength to the quantum system 
and to its measurement apparatus is application-dependent. When coupling a resonator to 
a qubit, either for a quantum memory or in a circuit quantum electrodynamics experiment, 
strong coupling to the qubit improves information transfer but also increases dephasing. 
When reading out a qubit, coupling the resonator strongly to its measurement apparatus 
increases the measurement bandwidth and signal but in addition increases dissipation [14]. 
Resonator designs therefore involve compromises between the competing needs for both 
strong and weak coupling [15, 16]. A resonator with a variable coupling would provide a 
significant improvement: If used to measure a qubit, the coupling to the measurement ap- 
paratus could be turned off except during resonator readout, when the coupling could be 
made large. When coupling two qubits through a resonator, the coupling could be turned 
on and off as needed [17, 18], yielding higher fidelity gates [4, 5]. 

Here we employ an externally- controlled variable inductance [19] to modulate the coupling 
of a resonator to a transmission line, creating the microwave equivalent of a Fabry-Perot 
cavity with a variable-transparency mirror. The resonator also has fixed coupling to a 
superconducting phase qubit. We demonstrate the time-controlled release of single-photon 
Fock and superposition states, thus generating a "flying qubit" [20-22]. We also perform 
timed capture and release of few-photon coherent states, and use the variable coupling to 
transmit and release photons with arbitrary waveforms [23, 24]. This new capability promises 
numerous applications in high-fidelity quantum computing and communication. 

The schematic in Fig. la displays a Fabry-Perot cavity, which represents the resonator, 
with a tunable transparency mirror to represent the variable coupler. A two-level atom 
plays the role of the qubit. In the actual experiment (Fig. lb, c and d), the resonator 
(r) is a quarter- wavelength (A/4) coplanar waveguide resonator, with one end coupled to 
a superconducting phase qubit {q) and the other end shorted to ground. Close to the 
grounded end (a distance ~ A/60 away), the resonator is connected to a variable coupler 
(c), which modulates the inductive coupling to a microwave transmission line. The vari- 
able coupler is controlled using a bias current, which flux- modulates the inductance of a 
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superconducting quantum interference device (SQUID) embedded in a mutual inductance 
circuit. The resonator-transmission line coupling Kc can vary from zero to a maximum emis- 
sion rate Kmax ~ 1/(5 ns), over a time scale of a few nanoseconds. The resonator frequency 
is fr — 6.57 GHz, and the phase qubit has a ground to excited state {\g) O |e)) tran- 
sition frequency tunable from ~ 6 to 7 GHz [7, 19, 25]. The qubit-resonator coupling is 
g/27r ~ 12 MHz, measured using swap spectroscopy [9]. The qubit-resonator interaction is 
controlled by tuning the qubit frequency, and is effectively turned off by setting the qubit 
frequency to its idle point, 400 MHz below the resonator frequency [7, 19, 25]. 

The resonator and variable coupler were characterized by measuring the decay of a one- 
photon Fock state stored in the resonator. When the resonator is weakly coupled to the 
transmission line, the photon decays due to internal resonator losses, while when strongly 
coupled, the photon is emitted into the transmission line. The pulse sequence is shown in 
Fig. 2a, where the qubit was excited from \g) to |e), and the excitation then swapped to the 
resonator, creating a one-photon Fock state [7] . The coupler bias current was then adjusted, 
and after a delay time r, the residual excitation swapped back to the qubit, and the qubit 
measured. 

Figure 2b displays the probability Pg of measuring the qubit in |e) as a function of delay 
T and the variable coupler current bias. Pg decays exponentially with time r, with the decay 
rate varying strongly with coupler bias. Two line cuts are shown in Fig. 2c, with exponential 
fits yielding the resonator lifetime Ti. For zero coupling, as determined from maximizing Ti 
with respect to coupler bias, we find the intrinsic Ti^j Ri 4.5 fis, in agreement with resonator 
loss measurements, while for coupling Kiarge the lifetime is reduced to Ti ^ 30 ns. The 
resonator inverse lifetime 1/Ti is the sum of the intrinsic decay rate 1/Pi,i and the coupler 
emission rate Kc, so Kc — 1/Pi — l/7i,i- The coupling dependence on current bias in Fig. 2b 
is in good agreement with calculations (Supplementary Information). 

We demonstrated dynamic control by changing the coupling during the delay period, as 
shown in Fig. 2d. We started with the coupling set to zero, and after a delay Tg switched 
the coupling to K\s.rge — 1/(30 ns). The reduction in the photon lifetime after the switch 
is clearly visible. The coupler switching speed was limited by the ~ 2 ns rise time of the 
coupler bias, roughly 2,000 times shorter than Ti^j. 

This measurement does not distinguish between incoherent decay and the expected phase- 
coherent release of the photon. We therefore also used heterodyne detection, with the 
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resonator "catching" and then "releasing" photons in coherent states. Figure 3a displays 
the pulse sequence: With the coupler set to an intermediate coupling Kc — 1/(356 ns), we 
excited the resonator with a 100 ns Gaussian pulse from the microwave source, with the pulse 
calibrated to trap (n) = 10 photons (see Supplementary Information). The coupling was 
then set to zero, trapping the photons for a storage time r^, then set back to Kc = 1/(356 ns), 
releasing the photons for heterodyne detection. 

Figure 3b shows the heterodyne-detected signal in the / (real) and Q (imaginary) quadra- 
tures in the time domain. During the Gaussian excitation pulse, the signal comprised the 
reflected component of the excitation. No signal was detected during the subsequent = 200 
ns storage time with the coupler turned off. A sharp onset followed by an exponentially- 
decaying envelope appeared when the coupler was turned back on, releasing the photons. 
The signal envelope has a decay time = 706 ns, in close agreement with the expected 
value 2/kc — 712 ns. The amphtude oscillations are from a 50 MHz mixer sideband signal, 
and the / and Q quadratures have a relative 7r/2 offset, as expected. As the output traces 
were averaged 10^ times, the presence of oscillations indicates that the output represents 
coherent photon release, with a fixed output phase relative to the input. 

Figure 3c displays the / quadratures using storage times r<j = 100 and 300 ns. These are 
identical during the state-generating Gaussian pulse, but during the release the oscillation 
phase depends on the storage time r^, scaling as ~ 0.8l7r(Ts/100 ns). This phase accumu- 
lation is as expected from the small tuning of the resonator frequency fr with coupler bias 
(see Supplementary Information), further demonstrating the coherence of the release. 

We also calculated the radiated energy f^''(P{t) + Q'^(t))dt, integrating the signal power 
from the photon release time tr to a cutoff tc = U + ST^. We find that the released energy 
for Tg — 300 ns is 4% lower than for = 100 ns, in agreement with the expected intrinsic 
resonator loss. 

These measurements confirm the phase-coherent capture and release of coherent states. 
To demonstrate that we can achieve the same control for non-classical states, we used the 
qubit to generate [7] the photon superposition state (|0) + \ l))/^/2 and measured the release 
signal after turning on the coupling (Fig. 4a). For an intermediate coupling Kc — 1/(320 ns) 
and a large coupling Kc — 1/(30 ns), the signal's exponential decay envelope has a time 
constant Td 625 ns and 69 ns, respectively, close to the expected 2/kc, verifying that the 
coupling determines the release rate. The integrated energy for intermediate coupling is 7% 
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lower than for large coupling, attributed to greater intrinsic loss from the slower release. 

We next tested the release and detection of the qubit-prepared superposition state 
cos(^^/2)|0) + e*'^sin(^^/2)|l), akin to previous work with static coupling [20, 21]. The re- 
leased photons were heterodyne-detected as a function of the Rabi angle 9 and the phase 
angle 0, with Fourier transforms of / and Q yielding the signal amplitude and phase. Fig- 
ure 4b shows the dependence of the signal amplitude on 9, with (j) — Q. The maximum 
amplitude is at 9 — 7r/2, corresponding to (|0) -|- |l))/-\/2. The amplitude goes to zero for 
the pure Fock states at ^ = and tt as expected, due to the loss of phase-coherence at the 
Bloch sphere poles. When releasing the state (|0) -|- e*'^|l))/v^ with 9 = 7r/2 and varying 0, 
the signal has constant amplitude and phase increasing linearly with 0. 

The on-demand, real-time gating of the coupler enables precise shaping of the photon 
release waveform. Figure 4c shows the tailored time-dependent release of the (|0) -|- |l))/-\/2 
photon state, modulating the coupling with a 200 ns Gaussian bias pulse with peak coupling 
Kp, followed by a 100 ns delay and then completing the release with Kc — 1/(320 ns). 
Figure 4c shows the / quadrature signal for Kp = 1/(320 ns), 1/(30 ns), 1/(10 ns) and 
1/(5 ns), with a Gaussian-like release waveform mimicking the coupler pulse. For the top 
three sub-panels, energy integrals show that 17.5%, 43.1% and 100% of the total stored 
energy is released during the pulse, with the remainder released after the 100 ns delay. For 
couplings Kp> 1/(10 ns), the release is completed during the Gaussian pulse. In contrast to 
fixed coupling, in which the waveform decays exponentially with time, the bottom sub-panel 
shows an exponential-like increase of the waveform, as needed for high-fidelity transfer of 
photonic information [4, 5]. 

We have demonstrated the phase-coherent, controlled capture and release of coherent and 
superposition photon states from a resonator, using a resonator-transmission line variable 
coupling. This powerful technique should allow long-range entanglement [9, 26-30], where 
the shaped release we display in the last experiment is a key ingredient for high-fidelity 
state transfer [4, 5]. This capability will further enable tunable coupling for resonator-based 
dispersive qubit readout, where time-domain control can minimize deleterious dephasing 
while maximizing measurement bandwidth and signal strength. 
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FIG. 1. Experimental setup, a, Schematic for cavity quantum electrodynamics with a variable 
transparency mirror, where Kc is the decay rate through the variable mirror and Kj the intrinsic 
photon decay rate. The cavity, two-level atom and variable transparency mirror represent the 
resonator, qubit and variable coupler, respectively, b. Experimental schematic: The left end 
of the resonator (r) is capacitively-coupled to a superconducting phase qubit (q) with coupling 
g/2Tr ~ 12 MHz. The resonator is 4.8 mm long with resonant frequency fr — 6.57 GHz and the 
qubit is tunable over the range from ~ 6 — 7 GHz. The resonator is connected to a transmission 
line through a variable coupler (c). c. The variable coupler is a transformer connected A/60 ~ 0.3 
mm from the grounded end of the resonator r, and comprises two fixed inductors Li = L2 = 480 
pH with a negative mutual inductance M = —138 pH, and a SQUID with tunable inductance 
Ls{^). The SQUID inductance is modulated using the flux from a current bias line through a 50 
pH mutual inductance, with a flux tuning of ~ ib<I>o/2 with ±25 fiA of bias current. Microwave 
excitation signals drive the resonator through the transmission line, and signals from the resonator 
are amplified and demodulated using a mixer driven by a local oscillator (LO). The demodulated I 
and Q signals oscillate at the LO sideband frequency, d. Micrograph of device (top), with details 
of phase qubit (bottom left) and variable coupler (bottom right). 
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Delay r (ps) Delay r (ns) 

FIG. 2. Characterization of the variable coupler with Fock state |1). a, Pulse sequence. Top: The 
qubit (red), de-tuned 400 MHz below the resonator, is excited to |e) by a vr pulse, then tuned into 
resonance with the resonator to perform an iSWAP, generating a one-photon Fock state. After 
a delay time r, a second pulse transfers the resonator state back to the qubit, and the qubit is 
measured using a projective single-shot measurement. Averaging 600 times yields the qubit excited 
state probability Pg- Middle: Coupler bias (green) used for b and c, starting with zero current bias 
followed by a rectangular pulse with variable amplitude setting the coupler strength during the 
delay period. Bottom: Coupler bias (green) for d, where coupling is switched from zero current to 
zero coupling (nzero) to Kiarge after a time Tg- b. Excited probability (color scale) versus delay r 
(vertical axis) and coupler bias in flux units (horizontal axis). Pe decays exponentially with r due 
to combined photon emission and intrinsic loss, with resonator lifetime Ti set by coupler bias, c. 
Vertical line-cuts of b display exponential decay of Pe, which gives resonator lifetime Ti. Resonator 
lifetime for zero coupling (Kzero) is intrinsic Ti^i = 4.5 /is, while for Kiarge it is reduced to Ti = 30 
ns. Inset shows Pe for short times for Kiarge- d, Pe versus delay r for coupling strength switching 
from zero to Kiarge- Decay rate switches from intrinsic lifetime (4.5 fis) to 30 ns, with transition 
taking ~ 2 ns. Blue, red, green and purple lines correspond to switching delays Ts of 200, 400, 600 
and 800 ns, respectively. 
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FIG. 3. Catch and release of photons in coherent states. Ah traces averaged 10^ times, a, Pulse 
sequence. Resonator (blue) is driven by on-resonance, 100 ns Gaussian excitation pulse (50 ns 
full-width at half-maximum (FWHM)), and coupler (green) tuned to Kc — 1/(356 ns). Gaussian 
pulse is calibrated so resonator catches a coherent state with (n) = 10 photons. Photons are stored 
for time Tg at zero coupling, then released with Kc = 1/(356 ns). b, Top sub-panel: Demodulated 
/ and Q quadrature signals for Tg = 200 ns, with 50 MHz sideband oscillations and a relative 7r/2 
phase shift (lines are guides to the eye). Signals include reflected part of excitation pulse, followed 
by release signal after delay r^, comprising a sudden onset with exponential decay. Bottom sub- 
panel: / on expanded scale, with sinusoidal fit (green line), comprising sideband oscillations with 
exponential decay envelope (dashed red line, time constant = 706 ns). The fixed phase with 10^ 
averages indicates phase coherence of photon release, c, / quadrature for trapping delays = 100 
ns and 300 ns, showing excitation pulse and a delayed photon release, with a delay-dependent 
phase shift. 
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FIG. 4. Storage, release, and heterodyne detection of non-classical photon states, a, Top: Su- 
perposition Fock state (|0) -|- |l))/\/2 is prepared with 7r/2 pulse to qubit, followed by iSWAP to 
resonator. Coupler is switched at t = from zero to either intermediate Kc — 1/(320 ns) or larger 
Kc — 1/(30 ns), showing faster decay for larger Kc- Middle and bottom: Heterodyne-detected I 
quadrature for both Kc values. Mixer LO sideband frequency is 110 MHz, with 10^ averages, b, 
Release of resonator superposition state cos(0/2)|O) -|- e*'''' sin(^/2)|l), with controlled Rabi angle 9 
and phase angle <j). After state preparation, variable coupling changed from zero to Kc = 1/(320 ns), 
and released photons heterodyne-detected. Fourier transforms of / and Q quadratures yield signal 
amplitude and phase (3.6 x 10^ averages). Top: Amplitude for fixed = has expected dependence 
on 9, with only superposition states generating a signal (Fock states have completely undefined 
phase). Bottom: Signal for fixed = 7r/2 as a function of 4>, with expected constant amplitude and 
phase scaling linearly with (j). c, First subpanel: Tailored, two-segment release pulse sequence for 
(|0) + |l))/\/2 superposition state, using a 200 ns Gaussian control pulse (100 ns FWHM, amplitude 
Kp) followed by a rectangular pulse with Kc = 1/(320 ns), with intervening 100 ns delay. Second 
to fifth subpanels: Pulse amplitude Kp set to 1/(320 ns), 1/(30 ns), 1/(10 ns) and 1/(5 ns). I 
quadrature amplitude depends on Kp, with different amounts of energy released during Gaussian 
pulse; two top sub-panels show remainder released during final rectangular pulse. 



12 



Supplementary Information for "Controlled catch and release of 

microwave photon states" 

Yi Yin/ Yu Chen/ Daniel Sank/ P. J. J. O'Malley/ T. C. White/ R. Barends/ J. 
Kelly/ Erik Lucero/ Matteo Mariantoni/'^ A. Megrant/ C. Neill/ A. Vainscnchcr/ 
J. Wenner/ Alexander N. Korotkov/ A. N. Cleland/'^ and John M. Martinis^' ^' * 

^Department of Physics, University of California, Santa Barbara, CA 93106, USA 
"^California NanoSystems Institute, University 

of California, Santa Barbara, CA 93106, USA 
^Department of Materials, University of California, Santa Barbara, CA 93106, USA 

"^Department of Electrical Engineering, 
University of California, Riverside CA 92521, USA 



* martinis@physics.ucsb.edu 



I. SAMPLE FABRICATION 



The device was fabricated using conventional photolithography and plasma etching on 
a multilayered structure. The resonator was made from a sputtered 150 nm-thick super- 
conducting aluminum base film on a sapphire substrate. The phase qubit and the super- 
conducting quantum interference device (SQUID) in the variable coupler were made using 
Al/AlOx/Al Josephson junctions. A low-loss dielectric, hydrogenated amorphous silicon, 
was used as the insulator in capacitors and wiring crossovers. 



II. EXPERIMENTAL SETUP AND CONTROL ELECTRONICS 

A detailed schematic of the resonator control and measurement system is shown in 
Fig. Sla. The microwave excitation signal for the resonator is generated by mixing a digitally- 
synthesized intermediate frequency (IF) signal with a microwave frequency local oscillator 
(LO). The I and Q quadratures of the IF signal were generated using a field-programmable 
gate array (FPGA) and a customized 2-channel 14-bit digital-to analog converter (DAC). 
The RF output of the mixer was sent into the cryostat and routed through a circulator to 
the variable coupler and resonator. Microwave signals from the resonator were routed by 
the circulator to a cryogenic amplifier G with 35 dB of gain and a noise temperature of 4 K. 
After further room-temperature amplification (60 dB of gain), the output signal was mixed 
down with the same local oscillator (LO) signal as the up-converter, generating the same IF 
(sideband) frequency. The / and Q quadratures were digitized at 500 megasamples/second 
using a 2-channel 8-bit analog-to-digital converter (ADC), and the signal passed to another 
FPGA for further processing. In the "oscilloscope mode", the digitized signals were sent 
directly to a computer without further processing, as shown in the time-trace data in Fig. 
3 and Fig. 4 in the main text. In the "demodulation mode" , the / and Q signals were mul- 
tiplied by sine and cosine waveforms at the sideband frequency and summed in real time. 
The continuous summation signal guarantees rapid fast Fourier transform (FFT) processing 
once the data acquisition is complete; this was used for the data in Fig. 4b in the main text. 

Qubit control used a system similar to that for the resonator. The microwave pulses for 
qubit control {x/y rotations) were generated by mixing a continuous microwave signal and 
a shaped quasi-d.c. waveform from a 2-channel 14-bit digital-to-analog converter (DAC). 
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Another quasi-d.c. pulse generator controls qubit z-axis rotations and measurement. 

Figure 1 in the main text shows the circuit for the experiment. The variable coupler is 
connected a distance d = 0.3 mm from the grounded end of the A/4 coplanar waveguide 
resonator (4.8 mm long). The coupler consists of a transformer with inductances Li, L2 
and a negative mutual inductance M. The positive mutual inductance from the dc SQUID 
is Ls — $o/47r/c| cos(7r$/$o)|) where Ic = 1.6 /lA is the critical current of the junction, $ 
the applied magnetic flux and $0 — h/2e the magnetic flux quantum. The characteristic 
impedances of the microwave resonator and the transmission line are Zr — 80 ft and Zq — 
50 Q respectively. 

III. THEORETICAL MODELING OF RESONATOR Ti 

Using an equivalent electrical circuit for the variable coupling experiment, shown in 
Fig. Sic, we calculate the expected resonator Ti due to coupling to the external 50 fl trans- 
mission line. We also calculate the small effect this coupling has on the resonance frequency 
of the resonator. Both calculations compare well with experiment. 

We replace the short portion (~ A/60) of the resonator between the coupler connection 
and the resonator ground with an effective inductance Lg, and the transformer and coupling 
circuit with an equivalent L[j L'2 and mutual inductance M', with L'2 connected to an infinite 
transmission line with characteristic impedance Zq. We calculate the resonator reflection 
and transmission amplitudes r and t. 

The effective inductances L\, L'^ and inductance M' are given by 

L; = Li + L„ L'2 ^L2 + L„ M' + L,. (2) 

The coupler is turned off when M' = 0. Note that because all the equivalent inductances 
include Lg, modulating M' by changing Lg modulates L[ and Lg as well: 

L; = (Li -M) + M', L'2 = (L2 -M) + M'. (3) 

To calculate the inductance Lg, which represents the small length of resonator from the 
coupling point to ground, we impose a voltage Ae^^^ in the resonator traveling from the 
open (qubit) end towards the coupler end, and approximate the reflected voltage as —Ae^^ 
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(ignoring the portion transmitted into the transmission hne) . The voltage at distance d from 
the grounded end isV ^ j^^iu^ti^^v^djc _^-u^dic^ ^ 2iAe^'^*sin(a;d/c), with c the phase velocity 
of light in the resonator. The current at this point is 7 = Ae'^\e''"^l'' jZr + e-^'^l'-jZ,) = 
2(A/Zr)e*'^* cos(a;(i/c), so the wave impedance is Z = V/I = iZrtsna^cjd/c) = icuLg, yielding 
the effective inductance 

Zr / ud\ Zr / 2nd\ 
(note this is evaluated at a; = 27r/r, the resonator frequency). 



The calculation of the transmission and reflection amplitudes t and r is similar to the 
derivation in [1]. Assume a voltage with amplitude A is incident on the coupler from the left 
side of the resonator, with reflected voltage rA and voltage transmitted into the transmission 
line tA. The voltage across L[ is V — {1 + r)A, while the voltage across L2 is denoted by 
X — tA. The currents flowing into L[ and L'^ are h — {1 — r)A/Zr — V/{iujLg) and 
I2 — —x/Zo, respectively. Using currents Ii and I2, we write equations for the voltage 
amplitudes x and V: 



X = iujM' 
[1 + r)^ = iuL'^ 



[l-r)A (1 + r)^ 



Z,- icuLe 
'l-r)A {l + r)A 



iujL'r, 



iuM' 



x 

X 



(5) 



From these equations, we can calculate the reflection amplitude r and transmission amplitude 
t = x/A (note that \t\^Zr/ZQ + jrp = 1): 



1 + r 



+ 



ZrZo{l+iujL'^/Zo) 



1 — r \ 

1 — a 



ZoLeil+iuiL'^/Zo) 



r = 



t =i 



1 + a 



2ujM' / I ia 
1 + a \Zr coLeJ l + iujL'2/Zo 



(6) 



In the limit ouLe <^ Zr and ouM' <^ Zq, which apply here, the reflection and transmission 
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amplitudes can be approximated as 



r — 1 + 2a — 1 + z- 



2uLeM' 1 
* Z,(L'i + Le) 1 + icuL^/Zo ■ ^ ^ 

The decay time of the resonator is obtained from the transmission amphtude [1] : 
_Q Zq _ 7rZrZo{L[ + L,)\l + u^L'i/Z^) 

where cu = 27r fr is the resonator frequency. 

The coupler bias dependence of the resonator Ti is extracted from the data in the main 
text and shown in Fig. Sic (blue dots). The predicted Ti from Eq. (8), using the actual 
circuit parameters, is also displayed in Fig. Sic (red line), in good agreement with the data. 
We note that the inductive coupling changes sign when the coupler strength sweeps through 
zero coupling [2, 3], verified by the expected tt phase change in a Wigner tomography mea- 
surement (see next section and Fig. S2d). 

The resonant frequency of the A/4 resonator is primarily determined by the resonator 
length and characteristic impedance, but is also affected by the variable coupler. The change 
in resonance frequency with coupler bias can be measured experimentally, and verified by 
the following calculation: Compared to the frequency fr at zero coupling (M' = 0), the 
resonance frequency shifts by A/, 

~ TT^Zr{Le + - M){Le + - M + M')' ^ ^ 

We compare the coupler bias dependence of fr, measured spectroscopically, with the fre- 
quency tuning from Eq. (9), in Fig. Sid. The frequency tunes over ~ 15 MHz, a very small 
fraction of the resonator frequency. Previous experiments have demonstrated resonator fre- 
quency tuning using Josephson junctions or SQUIDs embedded in a resonator [4, 5]. Here 
the frequency tuning is quite small, and is a by-product of the variable coupler located 
outside the resonator. 

We have measured the coupling dependence of the one-photon decay in Fig. 2b in the 
main text with a fast coupler bias. For comparison, a similar measurement was performed 



using the slow coupler bias, with data shown in Fig. Sle. The bias range is expanded to 
show the periodic response of the lifetime to the coupler flux bias. 

IV. CHARACTERIZATION OF RESONATOR STATES USING SWAP SPEC- 
TROSCOPY AND WIGNER TOMOGRAPHY 

We measured the single photon lifetime Ti to characterize the coupling strength, as 
discussed in the main text. We also used swap spectroscopy to perform an equivalent 
characterization, shown here with the coupler set to two representative coupling strengths. 
The pulse sequence is shown in Fig. S2a, starting with the system initialized in the ground 
state. The de-tuned qubit was excited by a tt pulse to |e) and then tuned close to resonance 
with the resonator, using a qubit tuning 2;-pulse with variable amplitude. The variable 
coupler was either left at zero coupling (^zero)) or switched to a coupling Kc — 1/(30 ns) 
immediately after tuning the qubit. In either case, the coupling strength was fixed for 
the full qubit-resonator interaction time r. The qubit excitation probability Pe was then 
measured using a triangular measurement pulse. 

The qubit excited state probability Pg is plotted versus the interaction time r and the 
qubit 2;-pulse amplitude in Fig. S2b and c, for weak and strong coupling, respectively. The 
chevron pattern due to the qubit-resonator photon swapping is evident in Fig. S2b, from 
which we calibrate the iSWAP [6] pulse amplitude and duration. In contrast, the response 
in Fig. S2c shows a rapid qubit-resonator relaxation, with energy strongly dissipated into the 
transmission line. We also note that the center of the chevron pattern in Fig. S2c shifts in 
comparison to Fig. S2b, due to the resonator frequency shift with coupler strength (Fig. Sid). 

The coupling strength changes sign when the coupler bias sweeps through the zero cou- 
pling point. An indirect phase-sensitive method, Wigner tomography [7], was used to detect 
this coupling sign change. The pulse sequence is shown in Fig. S2d. The resonator was 
prepared in the superposition (|0) + \ l))/\/2 state, and then driven by a variable-amplitude 
classical Gaussian microwave tomography pulse. The coupler was set to two different values 
during the tomography pulse, such that the coupling strength ±Kc had the same amplitude 
but opposite signs. The microwave tomography pulse, passing through the coupler, displaces 
and rotates the resonator state in the resonator phase space; the opposite coupling signs give 
opposite rotation directions to the resonator state for the same tomography pulse. The qubit 
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was then tuned on- resonance with the resonator for a variable time r, after which the qubit 
state was measured. Measurements of the qubit excitation probabihty Pe{T) were analyzed 
to yield the Fock state probability Pn{c^), where a is the complex amplitude and phase of 
the tomography pulse. The Wigner quasi-probability distribution VI/(q;) was calculated by 
evaluating the parity W{a) = ^(— 1)"P„(q;). 

The Wigner functions measured with the two signs of coupling strength ±Kc are shown 
in Fig. S2e. The Wigner functions clearly show a relative rotation angle of ~ tt. The density 
matrices of the resonator states can be calculated from the Wigner functions and projected 
onto the number basis pmn = ("^i|p|^), shown in the lower sub-panel of Fig. S2e. Here, 
we represent each element in the density matrix by an arrow, whose length and direction 
correspond to the magnitude and phase of pmn- The direction of the arrow for poi (pio) 
contains the relevant phase information, showing a ~ tt phase shift with a small phase error 
of 9.5°. 

V. CALIBRATION OF TRAPPED COHERENT STATE PHOTONS 

When a microwave Gaussian pulse with amplitude a and duration td is used to create a 
coherent photon state in the resonator, the state can be probed with a qubit through an on- 
resonance interaction, by measuring the qubit excited state probability Pe{T) as a function 
of the interaction time r. The photon state probability distribution Pn{ci) can be resolved 
in the Fock number basis \n) by decomposing Pe(T) into its discrete Fourier components 
fn = nfi, where /i = gf/vr is the vacuum Rabi frequency [7]. 

We used this measurement to calibrate the coherent state stored in the resonator for 
different coupling strengths and different microwave drive amplitudes a, with a fixed du- 
ration. In Fig. S3a we set the coupler strength to one of two values — 1/ (3000 ns) and 
Kc — 1/(210 ns), after which we excited the resonator with a variable amplitude Gaussian 
microwave pulse, and measured the qubit after a qubit-resonator interaction time r. The top 
sub-panel shows the pulse sequence, the middle sub-panci the qubit excitation probability 
Pg as a function of the microwave pulse amplitude and interaction time r for the smaller 
coupler strength, and the bottom sub-panel the same measurement for the larger coupler 
strength. A horizontal line cut (not displayed) shows a periodic but low amplitude oscilla- 
tion for small drive amplitude a, transforming to a clear ringing-collapse-revival pattern for 
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larger a. When the coupler is set to a small couphng (middle sub-panel of Fig. S3a), it is 
hard for the microwave source to excite the resonator but the resonator has a large Ti for 
trapped photons. When the coupler is instead set to a large coupling (bottom sub-panel of 
Fig. S3a), photons enter the resonator easily yielding a larger excitation amplitude, but the 
lifetime is shorter, illustrated by the rapid decay for larger r. 

To achieve both long photon lifetimes and low-power excitation, we instead set the coupler 
to a large value during the microwave drive pulse, then set the coupling to zero to trap 
the photons during the qubit measurement (Fig. S3b). A representative qubit- resonator 
interaction measurement for a coherent resonator state is shown in the middle sub-panel of 
Fig. S3b, with the coupler set to — 1/ (700 ns) during the microwave excitation pulse. We 
performed a series of measurements with varying coupling strengths, which were analyzed to 
give the photon distribution Pn(a). For a fixed microwave pulse amplitude (o; = 1.0 in the 
vertical axis of Fig. S3b), we calculated the average photon number {n{a)) — Ylim'^^nioi) . 
We display (n) as a function of coupler drive amplitude in the bottom sub-panel of Fig. S3b; 
this is the calibration method used for the experiment shown in Fig. 3 in the main text. 
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FIG. SI. Experimental setup and modeling of resonator lifetime Ti and resonance frequency fr- a, 
Full schematic for resonator control and measurement electronics. Qubit q is coupled with coupling 
strength 5 to A/4 resonator with characteristic impedance Z^- Other end of resonator is grounded, 
with variable coupler connected a distance A/60 from grounded end. Variable coupler comprises 
two inductances Li = L2 = 480 pH with a negative mutual inductance M = —138 pH, and a 
SQUID with tunable inductance -L<j(<I>). Current bias to the flux line to the SQUID controls the 
SQUID inductance and thus the variable coupler. The coupler is connected to a transmission line 
with characteristic impedance Zq, whose other end is connected through a circulator to a microwave 
excitation and measurement system (see text for details), b, Simplified schematic diagram of the 
variable coupler end of the resonator. The section of the A/4 resonator between the coupler and 
ground is approximated by an inductance Le, and the transformer is replaced by two effective 
inductances L'^, L'2 with mutual inductance M' . When a voltage signal with amplitude A travels 
from the left side of the resonator to the coupler, the wave is reflected from the coupler as rA and 
transmitted in the transmission line as tA. c. Blue dots are the experimental resonator lifetime 
Ti extracted from the data shown in the main text in Fig. 2b. Theoretical evaluation of Ti from 
Eq. (8) using the circuit design parameters is displayed as a red line, d. The resonator frequency fr 
as a function of coupler bias amplitude from spectroscopic measurements (blue dots) and compared 
with the theoretical prediction (Eq. 9). Arrows indicate coupling iLkc used for Wigner tomography 
in Fig. S2, with the sign of k denoting polarity of the inductive coupling, e. One-photon decay 
measurement, similar to Fig. 2b in the main text, but using a slow coupler bias. The periodic 
response of the lifetime to the coupler current bias is evident. 
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FIG. S2. Swap spectroscopy and Wigner tomography at different coupler biases, a, Pulse sequence 
for swap spectroscopy at two coupler settings. The qubit is excited to |e) with a vr pulse and the 
excitation swapped to the resonator. The coupler is either left at zero coupling when the qubit is 
tuned on- resonance with the resonator, or switched immediately to Kc = 1/(30 ns). After a qubit- 
resonator interaction time r, the qubit excited state probability Pg is measured, b, Measured 
qubit probability Pe for swap spectroscopy as a function of qubit z-pulse amplitude (detuning) 
and interaction time r, with the resonator coupling to transmission line set to zero. The chevron 
pattern representing qubit-resonator photon swapping is clearly visible, c, In contrast, when the 
coupler is set to strong coupling Kc = 1 / (30 ns) , swap spectroscopy shows a rapid energy dissipation 
and a slight resonant frequency shift, d. Pulse sequence for Wigner tomography. The box labeled 
"prep state" represents resonator preparation in the state (|0) + |l))/-v/2- The coupler is set to 
zero coupling during the entire sequence except when the microwave source drives the resonator 
for the tomographic analyzer pulse. During the tomographic pulse, the coupler is set to the same 
coupling strength but with opposite sign iLkc = ±1/(2000 ns). Following the tomographic pulse 
the qubit is used to measure the resonator state, e, Wigner functions lF(a) (upper sub-panels) 
for the (|0) + |l))/\/2 resonator state, plotted as a function of the microwave tomography complex 
amplitude a in photon number units (51 by 51 pixels). We calculate density matrices (lower sub- 
panels) from each Wigner function. The negative sign for the coupling strength introduces a tt 
phase shift between the tomography pulse-induced state rotations, with a small phase error of 9.5°. 
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FIG. S3. Calibration of coherent state generation for different coupler strengths, a, Top sub-panel: 
Pulse sequence to generate a resonator coherent state and then perform a qubit measurement. 
This measurement was performed for two different coupling strengths, Kc — 1/(3000 ns) and 
Kc — 1/(210 ns), with the coupler set to this value prior to the excitation pulse and left at this 
value during the qubit-resonator interaction. Data in bottom two panels show the qubit excited 
state probability Pg versus interaction time r and microwave drive amplitude. The Gaussian 
microwave pulse was 12 ns in duration (6 ns FWHM) for both panels a and b. b, Top sub-panel: 
Pulse sequence to generate a resonator coherent state with a coupling strength Kc — 1/(700 ns), 
with the coupling strength set to zero during the subsequent qubit-resonator interaction. Middle 
panel shows the qubit Pe as a function of interaction time r and microwave drive amplitude. A 
population analysis yields the average trapped photon number (n) for different coupling strengths 
during the microwave drive pulse, shown in the bottom sub-panel, for a microwave drive amplitude 
a = 1.0, the same as 1.0 in the vertical axis of the middle sub-panel. Coupler biases yielding zero 
average photon number are marked by arrows. 
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